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Abstract 

We discuss the interaction of photons with neutrinos including two 
lepton loops. The parity violation in the ■jv — > channel due to two 
\ lepton loops is substantially enhanced relative to the one lepton loop 

• contribution. However there is no corresponding enhancement in the 

parity conserving amplitude in either the direct or the cross channel 
6h! 77 -" 

> 

• r- 1 

X 

1 Introduction and Summary 

The photon-neutrino interaction is very feeble, involving only electrically 
neutral external particles. The cross sections are exceedingly tiny. Therefore 
this interaction can only be of astrophysical interest. Chiu and Morrison £Q 
proposed long ago that photon conversion to neutrinos may play a role in 
supernova cooling. Another possible application involves the propagation of 
light waves through a handed neutrino sea [2] which results in optical activ- 
ity (birefringence). There may perhaps be effects in the cooling of the early 
universe. For all these reasons there has been some small theoretical interest 
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in this interaction [3]. We consider here the interaction at energies which are 
small compared to the mass of the weak bosons. At these energies the weak 
interaction can be described by a Fermi type effective theory. To have z/7 
interaction one needs virtual leptons which couple both weakly and electro- 
magnetically, thus lepton loops in Feynman graphs. Contrary to intuition a 
two loop graph dominates the parity violation. In this section we summarize 
our results using only dimensional arguments. A more technical description 
is in the next sections. 

For supernova cooling the amplitude for the annihilation 77 — > vv is 
relevant, while for photon propagation in a neutrino sea the parity violation 
in the cross(scattering) channel 71/ — > is important. We discuss first the 
annihilation process where the two loop contribution turns out to be small. 
It is natural to expect higher order graphs to be smaller but as mentioned we 
find a violation of this rule here. This is reminiscent of a surprise in coherent 
photon scattering from atoms discussed by G.E.Brown and Woodward and 
later Peierls jl]. 

The earliest estimate for the annihilation amplitude used the Fermi theory 
of weak interactions |Tj , but Gell-Mann [5] noted that in V — A theory this 
amplitude vanishes for a point interaction of four fermions. This can be 
understood physically without any detailed computation. In Fermi Theory 
the annihilation amplitude is described by a triangle graph, with an electron 
running in the triangular loop (see Fig.l). 




Figure 1: 7^ scattering in the four-fermion Model. 

Two of the vertices of the triangle are electromagnetic (where the two 
photons couple) and the third vertex is weak, where the two neutrinos emerge. 
The amplitude has magnitude Gpa{pk), where p and k are the four- momenta 
of the neutrino and photon, and pk is the only non-vanishing relativistic 
invariant for the process (the Feynman amplitude is dimensionless which 
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brings the factor pk). Because the two neutrinos emerge from the same space- 
time point they can have no relative orbital angular momentum (s-wave) . 
With vector or axial vector coupling (V — A), the spins of the fermion-anti- 
fermion pair must be parallel. A theorem of Landau and Yang [S| forbids two 
free photons in a state with total angular momentum one. So the process 
is forbidden at order Gp in V — A theory. Gell-Mann also noted that this 
proof does not hold if there is a weak boson mediating the weak interaction, 
so that the two neutrinos are emitted at separate space-time points. 

With a weak boson to mediate the interaction, the (parity conserving) 
amplitude was estimated by Levine j7j , and his estimate remains valid in the 
standard model, as noted by a number of authors [3J IE] With the weak boson 
the loop changes from a triangle to a square, with the weak boson providing 
one of the sides (see Fig. 2). 




a) b) c) 

Figure 2: 71/ scattering in the Standard Model. 

This gives an additional factor of (pk/M^) to first order in the momenta 
p and k. This suppression factor vanishes when the boson mass goes to 
infinity, in agreement with Gell-Mann's theorem. The precise cross section 
for the annihilation was computed by Dicus and Repko [H] to order {Gfo) 2 : 

da/dz = [(G F a) 2 /32tt 3 ] [uj 6 /M^] [1 - z 4 ] , (1) 

where z is cosO and u the photon energy (in CoM system: pk = uo 2 (l — z)) 
and we have dropped logarithms. This is in agreement with the one loop 
Feynman amplitude described above, up to a constant factor: 

T"^ = [G pau 4 1 1 M^-][angular function][logarithms} . (2) 

The suppression factor (u 2 /My/) may be interpreted as a (kR) 2 factor as- 
sociated with a d-wave for the two neutrinos. Therefore the leading order 
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amplitude may be thought of as anomalously small compared to initial ex- 
pectations jS] • 

When we go to the next order in weak interactions there is a set of dia- 
grams in which the two neutrinos are emitted at a separation of the order of 
the Compton wavelength of the electron instead of the Compton wavelength 
of the weak boson. The diagram remains a square with two adjacent weak 
(order Gp) vertices, separated by an electron- neutrino loop. The two final 
neutrinos are emitted at the two weak vertices, and the two photons couple 
at the two remaining vertices (see Fig.3). 




Figure 3: The sample of two-loop contribution into 71/ scattering 

We can estimate the Feynman amplitude associated with the two-loop 
graph in lowest order using the same dimensional arguments as in the case 
of the one- loop result, changing the with m e in the denominator: 

T^ 2 ' = [G 2 F a(pk) s / m 2 ] = [G F au 6 /ml][angular function] (3) 

and therefore the ratio of the two-loop amplitude to the one loop is: 

T (2) /T (D = G F ^ { Ml/ml) . (4) 

The bracketed ratio of the two squared masses is very large ~ 10 10 so that 
the cross-section is dominated by the leading term only at low energies 
(u < m e ). However at intermediate photon energies u > m e our estimate 
of the two loop contribution is no longer valid since the ratio {pk)/m 2 is no 
longer small. To estimate what happens at these energies we take the limit 
where this ratio is large and then we expect m 2 e to be replaced by (pk) . Then 
the ratio of the two loop contribution to the one loop contribution becomes 
of order GpM^ which is a small number (of order 10 -2 ) independent of 
energy. Therefore there is no enhancement of the annihilation amplitude due 
to two loop diagrams. This is in accord with expectations. 
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We come now to the parity violating amplitude in the cross channel. 
This was first studied by Royer j2], using the triangle graph, and his result 
criticized by Stodolsky |2 on the basis of Gell-Mann's theorem. More recently 
estimates have been proposed using the standard model [3| and the one lepton 
loop graph, for both real and virtual photons. Here we show that these 
one loop estimates are all negligible when compared to the two lepton loop 
estimate-for real photons. 

In the cross channel the parity violating amplitude has to be odd under 
inversion of the coordinates. The precise form is discussed in section 2, and 
given below. As a result, parity violation appears dimensionally as a two 
loop effect even for a one loop amplitude. So it is natural to consider two 
loop terms together with one loop terms. 

The estimate of the parity violating Feynman amplitude in the standard 
model can be written in the form [3]: 

where C is a constant involving logarithms [3] and in the second line we 
evaluate the magnitude of the first expression in the forward direction for 
photons of energy u where pk = 2u 2 . This formula should apply for u < My/- 
With two lepton loops the corresponding parity violating Feynman am- 
plitude has the form(for u < m e ) 

T$ = 3G \^ r e, vaP e,(k)e:(k)[^) ~ 0.025G^ 6 /m*. (6) 

O.Z7T /TZg 

If we compare the two lepton loop to the one lepton loop result we find 

T$/Tff ~ (10-*)G F M 2 w [M 2 w /ml}. (7) 

This ratio is quite large and independent of the photon energy below 
m e ~ 500keV; its value is approximately 2 • 10 5 . Physically this reflects 
the fact that the two final neutrinos are separated in the two loop Feynman 
graph by an intermediate electron and a neutrino, and hence a distance of 
an electron Compton wavelength. Even for photon energies above lMeV the 
corresponding graph with a muon loop enhances the amplitude for photon- 
muon-neutrino scattering by a factor of 5. Should such an effect be observed 
it would indicate a violation of CPT and Lorentz invariance [5J. Attempts 
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to observe such an effect have been unsuccessful [TO]. Why is the situation 
different with the parity violation in the forward amplitude in the scattering 
channel? The answer is simple. This amplitude must be odd under parity 
and therefore be at least cubic in kp. So we need an extra factor of kp 
and this comes (in the one- loop case) with its own denominator Mfy, for 
dimensional reasons. Therefore the parity violating amplitude even in the 
one loop diagram is formally of order G 2 F and therefore on par with the two 
loop diagram. However it should be emphasized that even with this large 
enhancement factor the effect remains extremely small. 

For real photons k 2 = 0, the two loop estimates explained above, 

dominant. Therefore estimates of the rotatory power of the neutrino sea in 
the literature [3| have to be multiplied by a factor of 2 • 10 5 to make them 
right. In particular, using the Equation (31) of Abbasabadi and Repko [3], 
and applying this correction we get for the rotatory power <p/L the formula: 

<f>/L = 0M5G 2 F auj 2 T 2 {N„ - N p )/m 2 e , (8) 

which should only apply in the range discussed by Abbasabadi and Repko 
|3| ( where T v is the neutrino Fermi energy). 

For off-mass-shell photons the situation is different, since the Landau- 
Yang theorem does not apply to off-mass shell photons. Therefore even the 
lowest triangle graph is viable as was realized by Nieves et al [3] . At low pho- 
ton energy their one loop estimate remains dominant while at higher photon 
energies the two loop estimate dominates. The transition between Low and 
High photon energies u*, is given by the equation ui*E u > 30eV 2 , where E v is 
the neutrino energy in the rest frame of the sea. For a standard sea of energy 
of a few degrees Kelvin the transition occurs in the u ~ lOOifeV^region, but 
for an unconventional sea of Fermi Energy lOOeV^ jTT] the two loop estimate 
will dominate even for visible photons. For radio waves the one loop estimate 
dominates for virtual photons. 

2 Effective Lagrangian 

For a systematic study of neutrino-photon interaction it is useful to work 
with effective Lagrangians. 

The simplest example of the Effective Lagrangian is the four-fermion in- 
teractions of neutrinos v with electrons e. For fermions with momenta much 
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smaller than intermediate bosons mass one can integrate degrees of free- 
dom associated with W and Z and write the Effective Lagrangian only for 
fermionic degrees of freedom: 

G 

L eff = ^|(z>7a^)(er a e) , (9) 

where T a = gvla + QAlals- In the Standard Model (SM) one finds gy = 
| — 2 sin 2 9w, 9a — §■ For small momenta this effective Lagrangian is as good 
as the fundamental Lagrangian of the SM. 

Consider now the process v{p) + / y(k) — ► v(p) + / y(k) at momenta smaller 
than intermediate bosons mass or even smaller than the mass of electron m e . 
Interaction between neutrino and photons occurs only through the interac- 
tion of photons with charged virtual particles in loops. Loop diagrams are 
numerous and a bit complicated. For small momenta (pk/m,yy) <C 1 one can 
expand v(p) + j(k) — > u{p) + j(k) amplitudes in the power series in this 
small parameter. The lowest terms of this expansion can be represented as 
a matrix element of the appropriate operators i.e. of the particular terms of 
the Effective Lagrangian. 

Each term of the Effective Lagrangian has to be Lorenz-invariant com- 
bination of gauge-invariant electromagnetic field tensor F^ v and left-handed 
neutrino field v L = |(1 + 75)1/ and their derivatives. Effective Lagrangian 
has dimension four: [L] ~ [m] 4 . The operators may have higher dimensions 
D. To preserve correct dimension the coefficients in front of these operators 
should be proportional to appropriate power of (1/m), where m is the scale 
of mass walking inside the loops. The actual calculation of the diagrams 
gives numerical coefficient in front of the operator. 

This line of reasoning is very similar to the naive dimensional arguments 
given in the Introduction. The only advantage of effective Lagrangian is that 
within that more advanced approach we get more clear understanding of the 
structure of the operators (i.e. of the scattering amplitudes). 

2.1 P-even scattering amplitude 
2.1.1 One-loop approximation 

Consider effective Lagrangian for P-even z/7-scattering. The photon's part 
of amplitude should be even under parity, i.e. be the same for right-handed 
photons and left-handed photons. One has to construct appropriate Lorenz 
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invariant operators from the fields F^ v and v^. The combination of the fields 
that satisfies all these conditions looks as follows: 

e 4 

L eff liF^aF^V^adpil + + h.C. , (10) 

It has dimension D = 8. Matrix element of L e ff for forward scattering 
gives the amplitude 

T~^{pk)h{kye{k) , (11) 

If we identify parameter m in eq.(ll) with the largest mass in diagrams 
(i.e. with my/) we reproduce well known result [Zj up to the numerical 
constant. 

From this exercise it is absolutely clear that non-zero scattering ap- 
pears only in the second order in photon momenta, i.e. in order (pk) 2 . Thus 
one immediately concludes that any results of zero order in k, i.e. of order 
(Gfcx) ~ Oi 2 /mly (e.g. such as in ref. [IJG]), are erroneous. 

There is no way to violate Gell-Mann theorem within effective Lagrangians. 
To get the amplitude of the order Gf®. one needs operator with D = 6. By 
direct inspection one finds that there is no gauge invariant operator with 
D = 6. 



2.1.2 Two- loop approximation 

Consider now two-loop amplitudes with light particles (i.e. electrons and 
neutrino) in intermediate state between two external neutrino vertices. We 
expect that these diagrams are proportional to G 2 F a/m 2 . Thus to preserve 
correct dimension in effective Lagrangian we need operator of dimension 
D = 10. Appropriate effective Lagrangian is 

G 2 

L eff ~ -^-[F^F^Wv^dpd^l + 7 sH + h.c. (12) 

m e 

For this L e ff the scattering amplitude is of the third order in (pk) 

T~C-^{pk) 3 e{k)e*{k). (13) 

Thus for P-even scattering second order loops give correction of the order 
(GFM v/ )(pk/ml), i.e. small correction to the one-loop result. 
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2.2 Optical activity. P-odd scattering amplitude 

Now let us come back to P-odd effects in z/7 scattering and find appropriate 
operators in L e ff responsible for optical activity. The Lagrangian of dimen- 
sion D = 8 that depends on P- odd combinations of photon polarizations 
has the form 

Leff~—[F^Ftf][P7 a d{,(l + 'y B )v]+h.c. , (14) 
m w 

where F^ = ^e^^F^. 

The surprise is that this operator does not work in our case. On can check 
that the matrix element of F^F^p between photons with the same momenta 
and polarization ( forward scattering) is identically zero and this operator 
of D = 8 does not contribute into P-odd forward scattering. Thus P-odd 
effects are zero in (pk) 2 /m'y V order. The first nonzero effect is of the third 
order ~ (pk) 3 . Birefringence of neutrino sea is strongly suppressed! 

To find effect in the next order we have to look for the operators of higher 
dimension D — 10. One of these operators looks like follows 

L e ff -£[F tia (d y F fl i 3 ))[v'y ,df3d J (l + 75)1/] + h.c. (15) 

With this L e ff the forward scattering amplitude for a photon with mo- 
mentum k and for a neutrino with momentum p is equal to 

T = C{e A /M 2 ){pk/m 2 )h^ a ^{k)tl{k){p a kp/m 2 w ). (16) 

This amplitude has different contribution to left-handed and to right-handed 
photons scattering : Tll = —Trr. 

3 Actual calculations 

3.1 One- loop calculations. Real photons 

The actual calculation of the coefficient C has been done in one loop-approximation 
in [Hj with the results 

T = C(ey8n 2 s 2 )(pk/m 2 w ) 2 e^ a pelk)C(k)(p a kp/m 2 w ) , (17) 

where 

C = 4/3(ln(m^/m 2 )-ll/3), (18) 
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in the third reference in [3| and 



C = 4/3(ln(m^/m 2 )-8/3), 



(19) 



in the fourth one The reason for that discrepancy is unknown. Though 
numerically eqs.(18) and (19) differ by only a few per cent it would be inter- 
esting to understand whether there is a correct one-loop calculation. 

On the other hand we find that for a P-odd effect two-loop diagrams by 
many order of magnitude larger than one-loop contribution. Thus we can 
neglect any one-loop results. 

3.2 One-loop calculations. Off-shell photons 

The optical activity for off-shell photon was first considered by Mohanty, 
Nieves and Pal in They noticed that Gell-Mann prohibition theorem 
does not work for the off-shell photons. Thus one can expect that off-shell 
amplitude is of the first order in l/m 2 ^. Indeed 



where s = sinOw 

Eq. (20) differs from the original result of ref. jS] by factor 1/2. The 
reason is that triangle diagram was missing there. (This is the triangle dia- 
gram with Z boson and two photons. For real photons, k 2 = the triangle 
diagrams cancel each other in the SM since the SM is anomaly free. But for 
off-shell photons each triangle diagram gives a contribution proportional to 
k 2 /m 2 . The sum of all triangles is nonzero and the main contribution comes 
from the electron loop. This contribution has to be taken into account) 

3.3 Two-loop calculation 

The physical reason for the dominance of the two-loop diagrams under the 
one-loop in P-odd amplitudes is simple. To escape Gell-Mann's restriction 
one needs non-local interactions in order to include higher orbital momenta 
of the pair vv into annihilation process. In the one loop approximation 
the expectation value of the orbital momentum of the neutrino pair is ~ 
p/m^-The factor 1/mw measures the shortest separation of two neutrino 
during interaction (non-locality). 



T= (e 4 /8n 2 s 2 )e flua f 3 e[ t k)el(k)(p a k l3 /m 2 w )(k 2 /6m 2 ) , 



(20) 
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In the two-loop approximation the expectation value of neutrino pair 
orbital moment is ~ p/m e . The factor l/m e is due to (e~e + z/) in the inter- 
mediate states. 

Kinematically P-odd amplitudes are of the order of (pk) 3 (see section 
2.2). On dimensional ground we conclude that two-loop amplitudes are of 
the same order (pk) 3 (see secton 2.1.2). Thus one-loop P-odd amplitude has 
the same dependence on (pk) as the two-loop P-odd amplitude. Moving to 
the next order in electro- weak interaction we loose a small factor a/2ir , but 
win a great factor (rnly/m 2 ). The net effect is 

( T ( 2 )/T«) ~ (a/2ir)(m 2 w /m 2 e ) ~ 10 7 . (21) 
Actual calculation is rather lengthy. The result is 

= (13/27)(GV/647r 4 )(^ + g 2 A )(pk)\^e,(k)el(k)(p a kp/m 2 e ). (22) 
Thus the enhancement factor is 



T( 2 )/T« = (a/Q47rstn 2 e w )(m 2 w /m 2 e )(13/27C)(g 2 v + g 2 A ) ~ 10 5 , (23) 

where C is one-loop coefficient from eq.(8). We have lost two order of mag- 
nitude compared with naive estimate in eq.(10) mainly due to the large loga- 
rithmic coefficient C in one-loop amplitude. Still enhancement factor is very 
large ~ 10 5 ! 

This research was supported by NSERC-Canada and by RFBR grants 
00-15-96562. 
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